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The demonstration of isotropic polarization response from semiconductor quantum dots (QDs) is 
a crucial step towards the design of several optoelectronic technologies. Among many parameters 
that impact the degree of polarization (DOPf^i) of a QD system, the shape asymmetry is a critical 
factor. We perform multi-million-atom simulations to study the impact of the elliptical shapes on 
the electronic and polarization properties of single and vertically stacked InAs QDs. The comparison 
between a low aspect ratio (AR) and a high AR QD reveals that the electronic and the polarization 
properties strongly depend on the AR of the QD; the elongation of a tall QD allows tuning of the 
DOP[-^i over a much wider range. We then extend our analysis to an experimentally reported vertical 
stack of nine QDs (9-VSQDs) that has shown significant potential to achieve isotropic polarization 
properties. We analyse the contribution from the shape asymmetry in the experimentally measured 
large in-plane polarization anisotropy. Our analysis shows that the orientation of the base elongation 
controls the sign of the DOPr-^i; however the magnitude of the base elongation has only a very little 
impact on the magnitude of the DOPi^i . We further predict that the elliptical shape of the 9-VSQDs 
can only tune either DOP[no] or DOP[iio] f° r the isotropic response. Our model results, in agreement 
with the TEM findings, suggest that the experimentally grown 9-VSQDs has either a circular-base 
or a slightly [110] elongated base. Overall the detailed investigation of the DOP[^] as a function 
of the QD shape asymmetry provides a theoretical guidance for the continuing experimental efforts 
to achieve tailored polarization properties from the QD nano-structures for the design of optical 
devices. 

PACS numbers: 



I. INTRODUCTION 



Deployment of semiconductor quantum dots (QDs) in 
the active region of optical devices offers unique elec- 
tronic and optical properties which can be exploited to 
design several optoelectronic technologies ranging from 
lasers^ to semiconductor optical amplifiers (SOAsJpl or 
single photon sources^, where they have successfully over- 
come critical challenges such as extremely low thresh- 
old, high speed response, or entangled photon emission, 
respectively. However, in these applications, a critical 
design parameter is the polarization response of QDs, 
typically characterized in terms of eithe r de gree of po- 
larization [DOP = (TE-TM)/(TE+TM)PE1 or TM/TE 
ratid^, where TE-mode is measured along a direction in 
the plane of the QD, and TM-mode is measured along 
the growth [001] direction for the GaAs(OOl) QDs. Engi- 
neering of QD nanostrcutures to achieve isotropic polar- 
ization (DOP ~ 0) is critical for the implementation of 
several optoelectronic devices, for example semiconduc- 
tor optical amplifiers (SOA's). 

InAs QDs grown by the Stranski-Krastonov (SK) self- 
assembly growth process typically exhibit very poor po- 
larization response (DOP close to 1.0) due to the large 
compressive biaxial strain surrounding the flat shapes of 
the QDs. The strain induced splitting between the heavy 
hole (HH) and the light hole (LH) valence bands leads 
to a dominant HH character in the few top most valence 
band states, thus significantly suppressing the TM-mode. 
Therefore, previous studies of the single InAs QDs have 
reported very high values of the DOP, typically larger 



than 0.8.E" 

The polarization response of InAs QDs is influenced by 
several parameters such as crystal/atomic symmetry, QD 
shape, composition profile etc. The atomistic asymme- 
try of the underlying zincblende crystals implies that the 
[110] and [110] directions are inequivalent. This lowers 
the overall symmetry of a perfectly circular dome-shaped 
QD from Coo„ to C^. As a result, TE-mode in the plane 
of the QD does not remain symmetric and significant 
in-plane anisotropy may be observed even for an ideal 
circular-base InAs QDPSI Therefore, a single value of the 
DOP is not sufficient to characterize the polarizati on re- 
sponse of the QD systems. This, in the past studies^^, 
has lead us to define a direction-dependent value of the 
DOP, 
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where the direction, [n] = [110] or [110], associated with 
the DOPp^-] is same as the direction of the TEp^j-mode 
in the plane of the QD. 

The value of the DOPp^-j also strongly depends on 
the shape of the QDs, which is significantly affected by 
the growth dynamics of the self-assembly process dur- 
ing the growth of the capping layers and the post-growth 
annealing processes^. As a result, the shape of a SK 
self- assembled QD is far from being perfectly circular or 
square, as typically assumed in the past theoretical stud- 
ies of the polarization properties. Several experimental 
investigations have suggested that the actual shape of 
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the QDs significantly deviates from the ideal circular- 
base (for dome or lens) or square-base (for pyramid), and 
usually tends to elongate along the [HOp^Hial or along the 
[IlOpHU directions. 

Furthermore, recent advancements in the growth tech- 
niques have allowed to control the shape of QDs, leading 
to the fabrication of strongly elongated QD like nano- 
structures™ These offer an enhanced exciton oscillator 
strength and allow the realization of single exciton-single 
photon coupling to build the fundamental blocks for the 
solid state quantum information. 19 Such elongations of 
the QDs along the [n] = [110] or [110] direction can 
significantly alter the value of the DOP^j and may be 
exploited to achieve tailored polarization response for a 
desired operation. 

Brief overview of the past theoretical studies: 
Despite significant experimental evidence for the elonga- 
tion of the QD shapes and its prospective potential to 
tune the polarization properties, the impact of the base 
elongations on the value of the DOPp^j is only barely 
known. The previous theoretical investigations of the 
QD elongations are focused on the study and design of 
the fine structure splitting (FSS = energetic differen ce 
between the two bright excitons, ep; 10 ] - e^ioj^P^^"^ 
or the spin polarization^, with very little emphasis given 
to the study of the polarization properties ( compari son 
of the magnitudes of the TE and TM modes) .HHIIIZI 

Sheng et aZp3 applied an effective bond-orbital model 
and discussed the impact of the base elongations on the 
in-plane polarization anisotropy of the InGaAs QDs de- 
fined as: 
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and they concluded that the electron-electron interac- 
tions and the alloy intermixing have very little contri- 
butions in determining the polarization properties of the 
QDs. 

Mlinar et alW*, using an atomistic pseudo-potential 
model, focused on the InGaAs QDs and studied the im- 
pact of the (In,Ga)As alloy randomness on the polariza- 
tion properties. They concluded that the alloy compo- 
sition fluctuations can significantly change the in-plane 
polarization anisotropy and therefore the experimentally 
measured polarization anisotropy may not be considered 
as a reliable measure of the QD shape asymmetry. 

Schliwa et aZp2 based on their k.p calculations, stud- 
ied the impact of piezoelectricity on the QD optical prop- 
erties. Although they provided a detailed investigation 
of the electronic and optical properties of the dome and 
pyramidal shaped QDs as a function of their vertical as- 
pect ratio (height/base) , only one case of the pyramidal 
shaped QDs (series D in their paper) was investigated 
for the study of the base elongations (lateral aspect ra- 
tio). Furthermore, in their study of the lateral aspect 
ratios of the pyramidal QDs, they kept the overall vol- 
ume of the QDs unchanged by altering both the heights 



and the base diameters of the QDs. Since the QD ener- 
gies are strongly influenced by a change in their height 
parametei^, so the reported results did not isolate the 
impact of the QD base elongations. 

Nevertheless, a comprehensive quantitative analysis of 
the impact of the QD base elongations on the polariza- 
tion dependent room temperature ground state optical 
emissions still remain unavailable. This paper, therefore, 
aims to bridge this gap by providing a detailed study of 
the DOPp^j and Pol|| for the [110]- and [TlO] -elongated 
QDs. For the calculation of the polarization dependent 
optical modes (TE and TM), we take into account the 
highest five valence band states, instead of just a single 
top most v alenc e band state, in accordance with the re- 
cent studie d 5 * 10 ^ where it has been shown that the calcula- 
tion of the room temperature ground state optical spectra 
must involve multiple closely spaced valence band states 
to accurately model the in-plane polarizability and to 
avoid discrepancy between the theory and experiments. 
Our calculations show that despite tuning of the DOPn^i 
over a wide range, the elliptical shapes of the single QDs 
do not lead to an isotropic polarization. Therefore, we 
extend our study to multi-layer QD stacks. 

Isotropic polarization from multi-la yer stacks: 
Past theoretical^ and experimenta l 9 * 29 * 30 * studies have 
shown that the polarization response of the QDs can be 
drastically improved by growing large vertical stacks of 
QDs (VSQDs), consisting of many closely spaced QD lay- 
ers to exploit inter-dot strain and electronic couplings. 
A recent experimental study 9 -* has shown that a verti- 
cal stack of nine QDs (9- VSQDs) exhibits DOP[ 110 ] 
-0.fi: however the measured PL spectra showed a large 
anisotropy in the in-plane polarization modes: TEjy 10 ] 3> 
TE[ 110 i . our atomistic calculations^*, based on an assump- 
tion of ideal circular-base dome-shape for the QD layers, 
reported TErj 10 i 3> TE[ 110 ] in agreement with experi- 
mental PL spectra. We attributed this large in-plane po- 
larization anisotropy to a small increase in the TMfooii- 
mode (coming from an enhanced HH/LH intermixing) 
and a large decrease in the TE[ 110 ]-mode due to the [110]- 
oriented hole wave function confinements. 

A good qualitative agreement of our calculations with 
the experimental results assuming ideal circular-base for 
the QDs leads to a fundamental question that how much 
is the contribution from the QD shape asymmetry in 
the polarization anisotropy for the 9- VSQDs? This 
work, based on multi-million-atom calculations, provides 
the answer that the interfacial hole wave function con- 
finements have major contribution in the experimen- 
tally measured in-plane polarization anisotropy, which 
is counter-intuitive to the common notion where the 
shape asym metry is considered mainly response for such 
anisotropiespsEH 

Furthermore, as an [110]-elongation increases DOP[ 110 ] 
and reduces DOP[y 10 j , so we investigate the possibility to 
exploit it to balance the built-in in-plane anisotropy such 
that to achieve DOP [110] ~ and DOP [I10] ~ 0? This 
would lead to the design of the QD based SOAs inde- 
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pendent of the in-plane direction. Our study reveals an 
interesting property of the 9-VSQDs that its polarization 
response is very sensitive to the orientation of its elon- 
gation and both TE-modes can not be reduced below 
TM-mode. Therefore, either DOPjy 10 ] or DOP^io] can 
be tuned for an isotropic polarization, and not both of 
them simultaneously. 

Finally, the quantitative study of the elongation de- 
pendent DOPp^-j also helps us to determine the geometry 
of the 9-VSQDs. Our theoretical calculations accurately 
predict [110] elongation of the 9-VSQDs in the experi- 
ment, which is also consistent with the TEM images^. 

The remainder of the paper is organized in the follow- 
ing sections: section II defines QD geometry parameters 
and describes the three types of elliptical shapes that we 
study in this paper. Section III documents our method- 
ologies. Our results for two different AR single QDs are 
presented in sections IV-A and IV-B. Section IV-C is 
about the vertical stack of nine QD layers (9-VSQDs). 
The sections IV-A, IV-B, and IV-C are written as self- 
contained sections so that a reader interested in only one 
type of QDs may only require to read the correspond- 
ing section. Finally, we provide an overall summary and 
main conclusions of our results in section V. 



II. SIMULATED QUANTUM DOT SYSTEMS 

A. Geometry Parameters 

Figsjlja), (b), and (c) show three quantum dot geome- 
tries simulated in this study. The InAs quantum dots are 
embedded inside large GaAs buffers comprised of « 15 
million atoms (60x60x60 nm 3 ) for (a) and (b), and « 
25 million atoms (60x60x106 nm 3 ) for (c). The quan- 
tum dots are placed on top of 0.5 nm thick InAs wetting 
layers. 

We study three dome-shaped quantum dot systems: (i) 
a low aspect ratio (AR) InAs QD with 20 nm diameter 
(d) and 4.5 nm height (h), (AR = h/d = 0.225); (ii) a 
high AR InAs QD with d = 20 nm and h = 8.0 nm (AR 
= h/d = 0.40); (iii) a vertical stack comprised of nine QD 
layers (9-VSQDs) separated by 5 nm thick GaAs spacer 
layers, where each layer consists of an InAs QD with d 
= 20 nm and h = 4 nm. The geometrical paramete rs of 
this 9-VSQDs are taken from the recent experimental*^! 
and theoretical^ studies where it has shown great techno- 
logical relevance for achieving isotropic polarization re- 
sponse. 

In remainder of this paper, we label the single QD with 
AR = 0.225 as a "flat" QD and the single QD with AR 
= 0.40 as a "tall" QD. In a previous study^, it has been 
shown that the flat and tall QDs with an ideal circular- 
base exhibit drastically different electronic and polariza- 
tion properties. The hole wave functions tend to reside 
in the HH pockets for the tall QDs, when the AR > 0.25. 
This introduces a large anisotropy in the in-plane polar- 
ization (Pol||). Therefore, this paper analyses the impact 



of the base elongations for both types of the QDs. Fur- 
thermore, as the 9-VSQDs consists of strongly coupled 
QDs so it can essentially be considered as an extension 
of the single tall QD with a very large AR = 45/20 = 
2.25. Our calculations presented in the section (IV) con- 
firm that the two single QDs (flat and tall) exhibit dras- 
tically different polarization properties as a function of 
their base elongations, and the 9-VSQDs overall exhibit- 
ing many similar characteristics as that of the tall QD. 

Note that a typical SK growth of a large vertical QD 
stack generally results in an increase in the size of the 
upper layer QDs^, however, no such increase in the 
QD dimensions is reported in the experimental study^. 
Therefore, we keep the size of the QDs uniform for the 
9-VSQDs. 



B. Elongation along High Symmetry Axis 

In order to study the impact of the base elongations 
of the three QD systems described above along the 
high symmetry crystallographic directions ([110] and 
[110]) on their electronic and polarization properties, we 
consider three types of the elliptical shapes, defined by 
an elongation factor r\ — djuoj/dmQi which is a ratio 
of the QD lateral diameters along the [110] and [110] 
directions, for the above mentioned three QD systems: 

Type-I: As schematically shown in Fig. [ljd), the 
diameter d[no] (drj 10 i) is reduced by Ad for [110] ([110]) 
elongation, while keeping the other diameter drj 10 i 
(dfnoi) fixed at 20 nm. This will reduce the overall 
volume of the QD. 

Type-II: For this type of elongation, as schemati- 
cally shown in the Fig. [ije), we simultaneously change 
both d[no] and drj 10 ] diameters by equal amounts Ad 
(one diameter is increased by Ad and the other diameter 
is decreased by Ad). Once again, the volume of the QD 
decreases, however, at much slower rate compared to the 
Type-I case. 

Type-IIv: As a special case of the Type-II elonga- 
tion, we again change the QD diameters along the [110] 
and [110] directions similar to the Type-II case, but 
keep the overall QD volu me unch anged. Most of the 
previous theoretical studies^ 4 * 21 * 24 * ^ have only analysed 
this type of elongation, so it allows us to make a direct 
comparison with the existing results. For a circular-base 
dome-shape QD with d = 20 nm and h = 4.5 nm, the 
volume of the QD is (l/3)7rd 2 h cx d 2 oc 20 2 . Now to 
keep the volume unchanged, the other reasonable choice 
for the diameters djuo] and drj 10 i is 25 nm and 16 nm 
that results in d[no] xdrf 10 i = 400. Therefore for the 
Type-IIv elongation, we choose between 25 nm and 16 
nm for the diameters d[ 110 ] and dp; 10 i. 

As the volume of an ellipsoidal QD is proportional to 
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FIG. 1: Schematic of quantum dots are shown for theoreti- 
cal modeling, (a) A low aspect ratio (flat) dome-shaped InAs 
quantum dot with the base diameter and height of 20 nm 
and 4.5 nm, respectively, (b) A high aspect ratio (tall) dome- 
shaped InAs quantum dot with the base diameter and height 
of 20 nm and 8 nm, respectively, (c) A vertical stack consist- 
ing of nine InAs QDs (9-VSQDs), each with the base diameter 
and height of 20 nm and 4 nm, respectively. The geometry 
parameters are directly taken from the experimental study^. 

(d) Top view illustrating the base elongations for the Type-I 
elongation. The elliptical shape is formed by decreasing the 
diameter either along the [110] or along the [-110] direction. 

(e) Top view illustrating the base elongations for the Type- 
II and Type-IIv elongations. The elliptical shape is formed 
by simultaneously decreasing (increasing) the base diameter 
along the [110] direction and increasing (decreasing) the base 
diameter along the [-110] direction. For the Type-IIv elonga- 
tion, we select values for the d[no] and drr 10 i such that to keep 
the overall volume of the QD fixed, (f) Plots of the products 
of the QD diameters along the [110] and the [110] directions, 
d[no] XdrjiQi, as a function of the elongation factor (r\). For 
the fixed QD height, this product is directly proportional to 
the volume of the QD. 



the product of diameters along its major and minor axes 
(d[iio] and dpjQi), we plot this product as a function of 
the elongation factor 77 for the three types of elongations 
in Fig. [ljf), quantitatively showing the decrease in the 
QD volume for the three types of elongations. As it will 
become clear in the later sections, the large decrease in 
the QD volume (AV) for the Type-I elongation will sig- 
nificantly impact the electronic and the optical proper- 
ties, even dominating impact of the diameter changes. 
For the Type-II elongation, relatively much smaller de- 
crease in the volume will compete with the changes in 
the diameter and the net impact will be from Ad for the 
small values of rj, and from the AV for the large values 
of 77. 

It should be noted that in all of the three types of elon- 
gations, we keep the height of the QDs fixed. It has been 
showrPS that the electronic and optical properties are 
very sensitive to the height of the QDs, so by keeping it 
unchanged, we eliminate its contribution and only focus 
on the impact of the base elongations. We al so specify 
that the previous theoretical investigations^^!] of the QD 
elongations have used the term "lateral aspect ratio" for 
the ratio of the QD base diameters, which is equivalent 
to the elongation factor rj defined in this study. In order 
to avoid confusion, we use aspect ratio for the ratio of the 
QD height and base diameter; and use elongation factor 
for the ratio of the base diameters. Finally, by defini- 
tion, the elongation factor 77 is = 1.0, >1.0, and <1.0 for 
the circular, [110]-elongated, and [110]-elongated QDs, 
respectively. 



III. METHODOLOGIES 

The atomistic simulations are performed using NEMO 
3-D simulatoiSSS, which is based on strain energy 
minimization by using the valence force field (VFF) 
moc j e j37i38] anc j electronic structure calculations by solv- 
ing a twenty-band sp 3 d 5 s* tight binding HamiltoniarP^l. 
Both linear and quadratic piezoelectri c pote ntials are cal- 
culated by using the published recipe^^ and included 
in the Hamiltonian. The polarization dependent opti- 
cal transitions are computed from Fermi's golden rule by 
absolute values of the optical matr ix elements, summed 
over the spin degenerate states^ 2 '. The polarization de- 
pendent ground state optical transition intensity modes 
(TE[ 110 ], TErj 10 ], and TM| 001 ]) are computed as a cumu- 
lative sum of the optical transitions between the lowest 
conduction band state (el) and the highest five valence 
band states (hi, h2, h3, h4, and h5^. 

We want to emphasize here that all of the simulations 
are performed over very large GaAs buffers surrounding 
the QDs to properly accommodate the impact of long- 
range strain and piezoelectric potentials in the electronic 
structure and the optical transition calculations. For 
the strain relaxation, we use mixed boundary conditions: 
bottom fixed, periodic in the lateral directions, and the 
top free to relax. For the electronic structure calcula- 
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tions, we use closed boundary conditions. The dangling 
bonds at the surface atoms are passivated according to 
the published modeP^. 



IV. RESULTS AND DISCUSSIONS 

In the next two subsections, A and B, we present our 
results for the flat and tall QDs, respectively, and quan- 
titatively analyse the impact of the elliptical shapes on 
their electronic and polarization properties. The subse- 
quent subsection C presents results for the vertical stack 
of nine quantum dot layers (9-VSQDs). 

Factors that shifts the electron and hole en- 
ergies: Before we start our analysis of the three QD 
systems under investigation, we specify the factors that 
shift the electron and hole energies as the base diameters 
of the QDs are increased or decreased. We identify four 
major factors as follows: 

(i) Change in QD Volume (AV): Fig. [ijf) provides 
a quantitative estimate of the changes in the QD volume 
as a function of the elongation factor 77 for the three 
types of the elongations. Since the volume only decreases 
for the Type-I and Type-II elongations, so it will result 
in an increase of the electron energies and a decrease of 
the hole energies. 

(ii) Change in QD Diameter (Ad): The QD base 
elongations are based on the increase/decrease of the 
QD diameters along the [110] and [110] directions. 
While this decrease/increase in the diameters will have 
very little impact on the ground state electron energy 
el (due to its s-type symmetrical wave function), it will 
affect the electron and hole p-states energies due to 
their orientation along these directions. The increase 
(decrease) in the diameters will produce corresponding 
increase (decrease) in the electron energies and decrease 
(increase) in the hole energies, oriented along their 
directions. 

(iii) Strain: The strain directly modifies the band edges 
and thus impacts the electron and hole confinement 
energies. The electron energies are shifted by changes in 
the hydrostatic strain (e xx + e yy + e zz ) only, whereas the 
hole energies are affected by changes in both the hydro- 
static and the biaxial strain [e xx + e yy — 2e zz ). Simple 
analytical relations based on the deformation potential 
theory can be applied to estimate these changes^. 

(iv) Piezoelectric Potential: InAs/GaAs systems are 
strongly piezoelectric: the orientation and the magnitude 
of the piezoelectric potentials significantly impacts the 
orientation and the splitting of the electron and hole 
excited states. It should be noted that although the 
piezoelectric potentials do not directly shift the electron 
and hole p-state energies, but they determine the Ad 
induced changes by controlling the orientation of the 
p-states. 

Overall for the Type-I elongation, AV is large whereas 
Ad is small since we keep one diameter unchanged and 



only reduce the other diameter by Ad, so the impact 
of AV dominates. For the Type-II elongation, the QD 
volume only slightly decreases, so the impact of AV is 
small. However, we increase one diameter by Ad and 
decrease the other diameter by Ad, so the overall impact 
of Ad is much stronger. Finally, for the Type-IIv elonga- 
tion, AV=0, and Ad is roughly same as for the Type-II 
elongation. 



A. Flat Quantum Dot (AR=0.225) 

In this subsection, we study the impact of the elonga- 
tions on the electronic and polarization properties of a 
flat QD as shown by the schematic in Fig. [IJa) . The QD 
has a base diameter of 20 nm and a height of 4.5 nm (AR 
= 5/20 = 0.225). Such low AR QDs are more commonly 
obtained from the strain-driven SK self-assembly growth 
process and their electronic properties have been widely 
studied in the literature. 



1. Electronic properties of the flat QD 

Fig. [2] plots the lowest three conduction band energy 
levels (el, e2, e3) and the highest five valence band en- 
ergy levels (hi, h2, h3, h4, h5) for the three types of 
elongations: (a, d) Type-I, (b, e) Type-IIv, and (c, f) 
Type-II. The figures are plotted using the same scale to 
facilitate mutual comparison. In order to understand the 
shifts in the energies, we also plot wave functions and 
piezoelectric potentials for the few selected cases. Fig. [3] 
shows the top view of the wave function plots for the 
circular-base QD and for the few selected [110] and [110] 
elongations. The QD boundaries and dimensions are also 
marked in each case. 

Fig. [4] plots the total (linear+quadratic) piezoelectric 
potentials along the [110] (dotted lines) and [110] (solid 
lines) directions through the center of the QDs, «0.5 nm 
above their base. The quadruple nature of the poten- 
tials is clearl y eviden t which has been well established in 
the literatur e! 21 * 40 * 42 ! and is shown to strongly influence 
the orientation of the electron and hole p-states. In this 
case of the flat QD with AR=0.225, we find that the 
quadratic component of the piezoelectric potential does 
not fully cancel the linear component inside the QD re- 
gion in contrast to the previous k.p study^, where the 
quadratic and the linear components were found to fully 
cancel each other for the AR < 0.5. Since the piezoelec- 
tric potentials are a strong function of the QD shape and 
composition, so we find that their results cannot be gen- 
eralized to all types of QDs. Theoretical studies by Us- 
man et alW^ and Islam et ai.^ have also shown non-zero 
values of the net piezoelectric potentials for the similar 
QDs. 

Our calculations show that the piezoelectric potential 
plots have two peaks at the QD interfaces, one just out- 
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FIG. 2: (a, b, c) The lowest three conduction band energy levels (el, e2, and e3) are plotted as a function of the elongation 
factor (77) for the (a) TypeT, (b) Type-llv, and (c) Type-II elongations, (c, d, e) The highest five valence band energy levels 
(hi, h2, h3, h4, and h5) are plotted as a function of the QD elongation factor (77) for the (a) Type-I, (b) Type-Hv, and (c) 
Type-II elongations. The corresponding increase/decrease in the optical gap energy (E s ) is also specified in each case by using 
the vertical arrows. 
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FIG. 3: The top view of the wave function plots for the lowest three conduction band (el, e2, and e3) and the highest five 
valence band (hi, h2, h3, h4, and h5) states are shown for the circular-base QD and for the selected elongations of the QD. 
The intensity of the colors in the plots represent the magnitude of the wave functions, with the dark red color indicating the 
largest magnitude and the light blue color indicating the smallest magnitude. The boundaries of the QDs are also shown to 
guide the eye. 
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FIG. 4: The plots of the total piezo- 
electric potentials (linear+quadratic) are 
shown for (a) circular-base QD, (b, d, f) 
[110]-elongated QDs, and (c, e, g) [T10]- 
elongated QDs. In each case, the type 
and magnitude of elongation is specified. 
The solid red lines are plotted along the 
[110] direction through the center of the 
QD, 0.5 nm above its base. The dotted 
(broken) black lines are plotted along the 
[110] direction through the center of the 
QD, 0.5 nm above its base. The bound- 
aries of the QD region are also marked 
in each case by specifying the lengths of 
QD along the [110] and [TlO] directions, 
d[uo] and d [Tl0] . 



side the QD region and one just inside the QD region. 
The electron and hole wave functions are found to be 
more influenced by the piezoelectric potential peaks in- 
side the QD, which causes the lower electron p-state (e2) 
to align along the [110] direction and the ground hole 
state (hi) to slightly elongate along the [110] direction 
for the circular-base QD (see Fig. [3] for rj = 1.0). The 
alignment of the lower electron p-state (e2) along the 
[110] direction is in agree ment with the experimental re- 
ports for the similar QDsP 1 ^ 

Lowest three conduction band energies: From our 
atomistic relaxations, we find that the hydrostatic com- 
ponent (e xx + t yy + e zz ) of the strain remains unchanged 
for both the [110] and the [110] elongations. Since the 
electron energies are only affected by the hydrostatic 
strain, so the strain does not contribute in the shifts of 
the conduction band energies. Also shown in Fig. [4] that 
the piezoelectric potentials do not exhibit any significant 
change for the elongated QDs. The peaks outside the 
QD regions are only changed by 1 to 2 meV and the 
peaks inside the QD region are decreased by 2 to 4 meV. 
These relatively small changes in the potentials will not 
result in any noticeable effects on the e2 and e3 ener- 
gies. Therefore, we conclude that the electron energies 
are mainly affected by the changes in the diameters (Ad) 
and the volume (AV) of the flat QD, while the strain and 



the piezoelectric fields have only minor contributions as 
a function of 77. 

The lowest electron energy level (el) has s-type sym- 
metrical wave function and is mainly affected by the de- 
crease in the QD volume (AV). The QD diameter change 
(Ad) does not have any noticeable impact on its energy. 
In Fig.[2la), as the QD volume decreases for the [110] and 
the [110] elongations, a nearly symmetric increase in the 
el energy is calculated. When the volume of the QD is 
kept fixed as in Fig. [2jb) , el energy is almost unchanged 
confirming that Ad has only minor impact on el. Fi- 
nally, for the Type-II elongation (Fig. [2^c)) , AV is very 
small and a corresponding small increase in the el energy 
is observed. Fig. [3] shows that the wave function for the 
el state also retains its s-type symmetry, with only slight 
elongation along the direction in which the QD is being 
elongated. 

The excited conduction band states (e2 and e3) have 
p-type symmetry and thus gets strongly affected by all 
of the three types of elongations. The AV is once again 
a dominant factor, which pushes these energy levels to- 
wards higher values. This is evident from Fig.[2]ja) where 
both e2 and e3 increase in energy irrespective of the elon- 
gation direction. However, if AV=0 as in Fig. [2jb), the 
increase in the QD diameter reduces the energy of the 
state aligned along its direction and vice versa. 
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The type-II elongation is an interesting case where the 
two factors, AV and Ad, compete as the QD diameters 
along the [110] and [110] directions are changed by equal 
values. Since the lower p-state (e2) is always oriented 
along the major-axis of the elliptical shape, so any in- 
crease in the corresponding diameter tends to reduce its 
energy while the decrease in the QD volume pushes it 
towards the higher energies. For the small values of the 
elongation factor (ry = 0.67, 0.82, 1.22, and 1.5), e2 de- 
creases due to dominance of Ad induced shift. However, 
for the larger values of the elongation factor (77 = 0.54 
and 1.86), the AV induced upward shift overcomes the 
Ad induced downward shift and hence the energy of c2 
increases. 

The energy of the higher p-state, e3, always increases 
as a function of 77 because it is oriented along the shorter 
diameter direction and hence an increase in its energy is 
supported by both, AV and Ad. 

Electron p-state splitting: The energy difference 
between the p-states (Ae p = e3-e2) is an important pa- 
rameter of interest as it provides a measure of the confine- 
ment anisotropy between the [110] and [110] directions^ 
and is sometimes used to characterize the fine struc- 
ture splitting (FSS) 25 . We find that Ae p is always 
larger for the [110]-elongations as compared to the [110]- 
elongations for the same of 77. This is because for the 
circular-base QD (77 = 1.0), the cumulative effect of 
the underline zincblende crystal asymmetry, strain, and 
piezoelectricity results in Ae p ~ 2.4 meV and favours 
[110] direction for the e2 state (see Fig. [3] for 77 = 1.0). 
Therefore, any [110] elongation merely enhances this 
asymmetry, whereas a [110] elongation first needs to over- 
come this inherent ss2.4 meV splitting in order to flip the 
orientation of e2 state and hence results in overall lower 
values of Ae p . 

Separation between the lowest two electron en- 
ergies: Another parameter of interest for the laser de- 
sign is the difference between the lowest two conduction 
band energy levels (Ae2i = e2-el) which should be large 
to avoid undesirable occupancy of the excited states^. 
The largest reductions in Ae2i are calculated to be «1 
meV, «8 meV, and «10 meV for the Type-I, Type-IIv, 
and Type-II elongations, respectively. These small vari- 
ations in Ae2i suggest that the elongations of a flat QD 
does not deteriorate this parameter for the implementa- 
tion of laser operation. 

Highest five valence band energies: Figs. [2jd), 
(e), and (f) plot the highest five valence band energy 
levels (hi, h2, h3, h4, h5) for the three types of the 
elongations. The shifts in the valence band energy lev- 
els are more complicated to understand due to their in- 
termixed HH/LH characters and much stronger confine- 
ments within the QD region. The changes in their ener- 
gies exhibit drastic differences for the three types of the 
elongations. 

Due to the heavier effective mass and stronger confine- 
ment inside the QD region, the orientation of the hole 
wave functions are mainly determined by the piezoelec- 



tric potential peaks inside the QD region. In the case of 
the circular-base QD as shown in Fig. [3j hi and h2 are 
slightly elongated towards the [110]-direction due to the 
negative peaks of the piezoelectric potentials along this 
direction inside the QD region. 

For the [110]-elongations, the piezoelectric potential in- 
side the QD region slightly reduces and the decrease in 
d[Y 10 j diameter results in an enhanced impact of the larger 
negative peaks of the potential outside the QD along the 
[110] direction. These two factors favour forces the hole 
wave functions to align along the [110] direction as shown 
in Fig. [3j When the QD is elliptical with its major-axis 
along the [110] direction, the internal negative piezoelec- 
tric potential slightly increases and dominate to align the 
hole wave functions along the [110] direction. Overall, 
we find that the hole wave functions always tend to align 
along the minor-axis of the elliptical flat QD. 

The biaxial strain for the Type-I elongation slightly 
relaxes and along with large AV, it dominates the up- 
ward shift induced by Ad and pushes the hole energies 
towards lower values. This is clearly evident for the [110] 
elongations. However for the [110] elongations, since the 
hole wave functions exhibit stronger alignment along the 
[110] direction , the increase in energies coming from Ad 
is enhanced and therefore some of the hole energies can 
be observed shifting slightly upward in the Fig. [2](d) . 

For the Type-IIv elongation in Fig. [2je) , A V = and 
thus the impact of Ad, which is much larger than for 
the Type-I elongation, dominates. We also find that the 
biaxial strain relaxation is much weaker in this case, and 
thus its downward shift is also very small. As a result, 
the hole energies shift towards larger values. 

Finally, AV is very small for the Type-II elongation, 
and Ad is very large, so overall its impact dominates 
and shifts all the hole energies towards higher values. 
The impact of the biaxial strain relaxation is again very 
small downward shift. 

As a summary, for the flat QD under elongations, the 
shifts in the electron and hole energies are mainly gov- 
erned by AV and Ad, whereas the direct impact of the 
strain remains negligibly small. 

Optical gap energy, E g : The optical gap energy 
(E 9 = el-hl) increases for both, the [110] and the [110], 
Type-I elongations, thus blue shifting the ground state 
optical wavelength mainly due to a large decrease in the 
QD volume. However, if the QD volume is fixed or only 
slightly decreased as for the Type-IIv and Type-II elon- 
gations, respectively, E ff decreases as a function of the 
elongation and hence results in a red shift of the ground 
state optical transition wavelength. 

2. Polarization properties of the flat QD 

Figs. [5](a), (b), and (c) compare polarization depen- 
dent TE and TM modes for the three QD elongations 
under study. The elliptical shape of the QD, irrespec- 
tive of the orientation of its major-axis, tends to increase 
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FIG. 5: The polarization dependent optical transition modes TE[no], TE[j 10 i, and TMpoi] are plotted as a function of the (a) 
Type-I, (b) Type-IIv, and (c) Type-II QD elongations. The figures are plotted using the same scales to facilitate easy mutual 
comparison, (d) The plots of the in-plane polarization anisotropy (Pol| | ) as defined by Eq. [5] are shown for the three types 
of elongations, exhibiting an inverse quadratic dependence on 77, consistent with the Sheng et a/. 26 . Four cases are marked 
by using green ovals indicating that the two different types of elongations, with roughly similar values of 77, exhibiting similar 
values of the P0I11. 



the TE-mode along its major-axis and decreases the TE- 
mode along its minor-axis. This can be understood as 
follows: the few top most valence band states have domi- 
nant heavy hole (HH) character due to the strain induced 
large splitting between the HH and LH bands. These 
heavy hole states are mainly comprised of \X) and |F) 
symmetry wave functions, where X and Y are selected 
along the high symmetry [110] and [110] directions, re- 
spectively. The lowest electron state (el) is mainly sym- 
metric IS") type wave function (as a good approximation). 
The elongation of QD along, for example, the A-direction 
will have negligible impact on the IS") type wave function, 
but it will increase (decrease) \X) (\Y)) component of 
the valence band states. Therefore, the TE^ cx |(X|S')| 2 
component of the electron-holes transition will increase 
and the TEy oc | | 2 component will decrease as evi- 
dent from Figs.[5ja), (b), and (c). 

The analysis of the calculated TMrnoi] component re- 
veals that it also increases for the elliptical QDs. The 
dominant contribution in the ground state optical in- 
tensity for the flat-shaped QDs comes from the high- 
est valence band state (hi), with the lower valence band 
states (h2, h3, h4, and h5) only adding weak transition 
strengths^. The strength of the TMrgoi] component is di- 
rectly related to the LH mixing in the valence band states 
and its magnitude is proportional to the electron-hole 
wave function spatial overlap along the growth ([001]) 
direction, given by elhi = |(?/'ei|^|V'fti)l- Here ipei is the 
ground electron state, ijihi is the ith valence band state 
where i € { 1, 2, 3, 4, 5}, and z is along the [001]- 
direction. 

Table [I] provides values of the ratios of the LH com- 
ponent in the top five valence band states for the few 
selected elliptical shapes with respect to the circular- 
base shape. For each case, the corresponding ratio of 
the spatial overlap (elhi) is also provided with in the 
(). Our calculations show an increasing LH mixing in 



the valence band states for the elliptical QDs with re- 
spect to the circular-base QD, in particular for the hi 
state which gives dominant el-hl transition. The spatial 
overlap between electron and hole wave funcitons also in- 
creases and therefore, a net increase in the TM[o i] mode 
is calculated as a function of 77. This characteristic of the 
pure InAs QDs is in contrast to the In 5 Gao.sAs ordered 
and disordered QDs reported by Singh et al^, where 
they have shown that the LH character of the hi valence 
band state remains unchanged as a function of the QD 
elongation. Based on these results, we predict that the 
elliptical shape has a stronger impact on the polariza- 
tion response of the pure InAs QDs, as compared to the 
alloyed InGaAs QDs. 

In-plane polarization anisotropy: Fig. [5jd) plots 
the in-plane polarization anisotropy (Pol| | ) as defined by 
Eq.[2]for the three types of the elongations. We find that 
our results for Pol|| are overall in agreement with the find- 
ings of Sheng et alW^, that when the height of the QD is 
kept fixed, similar elongations (77) results in nearly similar 
values of the in-plane anisotropy. We highlight four such 
cases in Fig. |5jd) using green dotted circles where nearly 
same values of Pol| | are calculated for the following sets 
of values of 77: (1) 0.64 in Type-IIv and 0.67 in Type-II, 
(2) 0.82 in Type-II and 0.8 in Type-I, (3) 1.22 in Type-II 
and 1.25 in Type-I, and (4) 1.565 in Type-IIv and 1.50 
in Type-II. Therefore, we conclude that the in-plane po- 
larization anisotropy (Pol||) only depends on the value 
of 77, irrespective of the type of the elongation provided 
that the height of the QD is kept fixed. We also find that 
our calculated values of Pol| | for all of the three type of 
the elongations roughly follow an inverse quadratic de- 
pendence on the elongation factor (77) , in agreement with 
the quadratic dependence on the lateral aspect ratio (/3) 
reported by Sheng et a?P3 as by definition 77 = /3 _1 . 

Tuning of DOP^: Figs. |6ja,) , (b), and (c) inves- 
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TABLE I: The ratio of the LH components in the highest five valence band states (hi, h2, h3, h4, h5) of the elliptical-shaped 
flat QD with respect to the LH components of the corresponding valence band states of the circular-base (r/) QD for a few 
selected values of the elongation factor, 77. For each case, the corresponding ratio of the electron-hole wave function spatial 
overlap along the 2 = [001] direction given by elhi = |(V , ei|^|'0hi}| is also provided. 
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FIG. 6: The plots of the DOP[nn] and DOP[j 10 ] are shown as a function of the (a) Type-I, (b) Type-IIv, and (c) Type-II QD 
elongations. The values of the DOP decrease irrespective of the direction of the QD elongation. 



tigate the changes in the DOPp^j for the three types of 
the elliptical shapes. Significant changes in the values of 
the DOPp^-j are observed as a function of 77. More inter- 
estingly, it is found that both DOP[n ] and DOP[y 10 j 
decrease irrespective of the type and direction of the 
elongation. This implies that the elliptical-shape, in 
general, improves the polarization response compared to 
the circular-base for the flat QD. The largest decrease 
(«23%) in the value of the DOP [U0 ] (from 97% to 74.5%) 
is calculated for the [110] elongation (77 = 0.54). 

B. Tall Quantum Dot (AR=0.40) 

In this subsection, we study the impact of elliptical 
shapes on the electronic and polarization properties of a 
tall InAs QD, having the same base diameter (20 nm) 
as of the fiat QD of the previous subsection, but with a 
height of 8 nm (AR = h/d = 8/20 = 0.40). Such high 
AR QDs are obtained using special growth conditions 
(very slow growth rate and high temperatures)^, or are 
typically found in the optically active upper layers of the 
weakly coupled bilayer QD stacks 5 , where the presence of 
strain from the lower QD layers results in the larger size 



of the upper layer QDs. To our knowledge, no detailed 
theoretical investigation of the polarization properties of 
such tall elliptical dome-shaped QDs is available in the 
literature, as the previous theoretical studies have only 
focused on the flat shapes of the QDs with low ARs: AR 
= 2/20 = 0.]P1 AR « 0.1/fSl, AR = 3.5/20 = 0.175^1, 
AR = 2/25 = 0.08 27 , and AR = 4.5/28.8 = 0.HPI 

Schliwa et alW^ applied k.p theory to study the elec- 
tronic and optical properties of the QDs by varying their 
AR from 0.17 to 0.5. However in order to vary the AR, 
they choose to keep the QD volume constant by simulta- 
neously changing both the base diameter and the height 
of the QD. Therefore their results do not isolate the im- 
pact of only varying the lateral aspect ratio. Their study 
of the base elongations which is focused at the pyramidal- 
shaped QDs shows that the linear and quadratic piezo- 
electric potentials fully cancel each other for all values 
of the AR < 0.5, whereas our atomistic simulations pre- 
sented in the previous subsection have already shown a 
non-zero net piezoelectric potential in the interior of a 
dome-shaped QD with the AR = 0.225. We therefore 
conclude that their result can not be generalized for all 
types of the QDs. In order to extend our study of the pre- 
vious subsection, here we investigate a tall QD by just 
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increasing the height of the QD from 4.5 nm to 8 nm, 
keeping its base diameter same (20 nm) as for the fiat 
QD. This allows us to make a direct comparison between 
the results for the flat and the tall QDs. Once again, for 
all types of the elongations, the height of the QD is kept 
constant at 8 nm. 



1. Electronic properties of the tall QD 

Fig. [7] plots the lowest three conduction band energies 
(el, e2, e3) and the highest five valence band energies 
(hi, h2, h3, h4, h5) as a function of the elongation factor 
(77) for all of the three types of elongations: (a, d) Type-I, 
(b, e) Type-IIv, and (c, f) Type-II. In order to develop 
an understanding of these energy shifts, we also plot the 
wave functions and the piezoelectric potentials for a few 
selected values of 77 in Figs. [8] and [9j respectively. 

The piezoelectric potential plots in Fig. [9] show famil- 
iar quadrupole symmetry, but exhibit three noticeable 
differences when compared to the plots for the flat QD 
(Fig.§: 

(i) The interior of the QD is nearly field free, with only 

one peak at the QD interfaces. This is in agree- 
ment with Schliwa et aiP^ indicating a cancellation 
between the linear and the quadratic components 
inside the QD region. 

(ii) A much larger magnitude of the fields (nearly twice) . 

(iii) A flip in the sign of the fields. The potential peaks 
are positive along the [110] direction and negative 
along the [110] direction. 

It is interesting to note here that although the piezo- 
electric potential profiles are drastically different for the 
flat and the tall QDs, but the electron and hole states 
for the circular-base case (77 = 1.0) align in the same di- 
rection for both cases: e2 along the [110] direction and 
all hole states along the [110] direction, as shown in the 
Fig. [8) Thus, an experimental measurement on a circular- 
base dome-shaped QD would not indicate any sign dif- 
ference for the e[ 110 ]-e[Y 10 j as a function of QD AR, de- 
spite the underlying physical details are quite different 
and require atomistic modeling with realistic simulation 
domains and physical parameters. 

We also find that the wave functions of the hole states 
for the tall QD are localized close to the QD interfaces 
due to the presence of the heavy hole (HH) pockets in 
the valence band edges. In our previous studjffSl we have 
shown that such interfacial localization of the hole wave 
functions for a pure InAs QD starts when the AR is in- 
creased above 0.25. Similar results were presented by 
Narvaez et apHfor the InAs QDs using pseudo-potential 
calculations, where they varied the QD AR from (5/25.2) 
0.198 to (7.5/25.2) 0.298 and showed that the hole states 
tend to confine in the HH pockets at the QD interfaces 
for the tall QDs. 



Lowest three conduction band states: As for the 

flat QD of the previous subsection, the hydrostatic strain 
component does not change as a function of 77, thus the 
strain contribution in the conduction band energy shifts 
is negligible because they are only affected by the hydro- 
static strain. The piezoelectric fields also only slightly 
change when the QD is elongated (Fig. [9]) , so the contri- 
bution from the changes in the piezoelectric fields in the 
shifts of the electron energies is minor, with dominating 
contributions coming from AV and Ad. 

The lowest conduction band state el has s-type sym- 
metry and the impact of Ad on its energy is negligible. 
The Type-I elongation considerably reduces the QD vol- 
ume (see Fig. [TJf ) ) , and therefore el energy increases. 
When the QD volume is unchanged as in the case of the 
Type-IIv elongation or is only slightly decreased as in 
the case of the Type-II elongation, el shows a very small 
change. The wave function plots for the el state indicate 
s-type symmetry with only slight elongation, mainly for 
the [110] oriented Type-IIv and Type-II elongations. 

The excited electron states, e2 and e3, are separated 
by «8.8 meV for the circular-base QD (77=1.0), which 
is around four times larger than the «2.4 meV splitting 
for the flat QD. This larger splitting is mainly due to 
the larger (nearly twice) magnitude of the piezoelectric 
potentials as evident from the comparison of Figs. [4] and 

M 

For the Type-I elongation (Fig. [7][a)), the large re- 
duction in the QD volume shifts both, e2 and e3, to- 
wards the higher energies. The energy difference Ae p 
= e3-e2 increases for the [Tl0]-elongation and decreases 
for the [110]-elongation, which is in agreement with 
what Schliwa et alW^ also calculated for the square-base 
pyramidal-shaped QDs. They reported electron p-state 
degeneracy towards the [110] base elongation. 

When the overall QD volume is fixed in the [110] Type- 
IIv elongation, the Ad induced shift lowers the energy of 
the [110]-oriented e2 and increases the energy of the [110]- 
oriented e3, as shown in Fig. [7^b). Since e2 and e3 go 
through a flip of their orientations for the [110] Type-IIv 
elongation, so only a very small change in their energies 
is observed (from 8.8 meV to 9.3 meV). 

Finally, the changes in the energies of e2 and e3 for the 
Type-II elongations (Fig. [7j[c)) follow the similar trends 
as previously calculated for the flat QD (Fig. [2]jc)). The 
two competing factors (AV and Ad) contribute to the 
energy shifts, where Ad being dominant for the small 
elongations (77 = 0.67, 0.82, 1.22, and 1.5) and AV taking 
charge for the large elongations (77 = 0.54 and 1.86). 

Separation between the lowest two electron en- 
ergies: The difference between the lowest two conduc- 
tion band energy levels (Ae2i = e2-el) which is im- 
portant for the laser design robustness increases for the 
Type-I elongations, with largest increase being «8 meV 
calculated for the [110]-elongation. For the Type-II and 
Type-IIv elongations, only minor reductions in the Ae2i 
are calculated. Therefore, we extend our conclusion of 
the previous subsection that the elongations of both, the 
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FIG. 7: (a, b, c) The lowest three conduction band energy levels (el, e2, and e3) are plotted as a function of the elongation 
factor (77) for the (a) Type-I, (b) Type-IIv, and (c) Type-II elongations, (c, d, e) The highest five valence band energy levels 
(hi, h2, h3, h4, and h5) are plotted as a function of the QD elongation factor (77) for the (a) Type-I, (b) Type-IIv, and (c) 
Type-II elongations. The corresponding increase/decrease in the optical gap energy (E 9 ) is also specified in each case by using 
the vertical arrows. 
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FIG. 8: The top view of the wave function plots for the lowest three conduction band (el, e2, and e3) and the highest five 
valence band (hi, h2, h3, h4, and h5) states are shown for the circular-base of the QD and for the selected elongations of the 
QD. The intensity of the colors in the plots represent the magnitude of the wave functions, with the dark red color indicating 
the largest magnitude and the light blue color indicating the smallest magnitude. The boundaries of the QDs are also shown 
to guide the eye. 
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FIG. 9: The plots of the total (lin- 
ear+quadratic) piezoelectric potentials 
are shown for (a) circular-base QD, (b, 
d, f) [110]-elongated QDs, and (c, e, g) 
[110]-elongated QDs. In each case, the 
type and the magnitude of the elonga- 
tion is specified. The solid red lines are 
plotted along the [110] direction through 
the center of the QD, 0.5 nm above its 
base. The dotted (broken) black lines are 
plotted along the [110] direction through 
the center of the QD, 0.5nm above its 
base. The boundary of the QD region 
is also marked in each case by specifying 
the lengths of QD along the [110] and 



[110] directions, d[ 110 j and d| 



no]- 



flat and the tall QDs, do not deteriorate Ae2i for the 
implementation of laser operation. 

Highest five valence band energies: While the 
conduction band energy level shifts for the tall QD quite 
resemble with the corresponding shifts for the flat QD, 
the hole energy level shifts show significant contrasts. 
The main reason for this different behavior is the dis- 
similarity of the net piezoelectric potential profiles for 
the two types of the QDs as evident from the comparison 
of Figs. [4] and [9j which leads to very different confine- 
ments of the hole wave functions (see Figs. [3] and [8]). 
The large negative peaks of the potentials at the QD 
interfaces along the [110] direction and nearly field free 
interior of the tall QD result in the hole wave functions 
being oriented along the [110] direction irrespective of 
the type and the direction of the elongations as shown 
in Fig. [8] Furthermore, the HH pockets at the QD in- 
terfaces results in the hole wave function confinements at 
the QD interfaces. 

For the Type-I elongation (Fig. [7][d)), both a large de- 
crease in the QD volume (AV) and a small relaxation of 
the biaxial strain pushes the hole energy levels towards 
the lower energies. When the QD is elliptical along the 
[110] direction, the hole wave functions also being ori- 
ented along this direction experience a downward shift in 
their energies. Thus all of the three factors add to result 
in a stronger reduction in the hole energy levels. In the 
case of the [110] elongation, d[ 110 ] is reduced which pushes 
the hole energy levels upward. Although the cumulative 



downward shift from AV and biaxial strain relaxation 
overall remains dominant, however for some values of 77, 
Ad induced upward shift is evident. 

As AV=0 for the Type-IIv elongation and the biax- 
ial strain relaxation is much weaker as compared to the 
Type-I case, so only Ad induced shifts are observed in 
Fig. Re): the increase in d[ 110 ] decreases the hole ener- 
gies (for the [110] elongation) and the decrease in dr 110 i 
increases the hole energies (for the [110] elongation). 

Finally, for the Type-II elongation in Fig. [7]T), AV is 
quite small. The biaxial strain relaxation again causes 
a small downward shift in the hole energies. However, 
the dominant shift comes from the Ad. Therefore, the 
hole energies move towards the lower values for the [110] 
elongation and shifts towards the higher values for the 
[110] elongation. 

Optical gap energy, E g : The changes in the optical 
gap energy, E g — el - hi, are also marked with the help 
of vertical arrows (using red color) for the three types of 
the elongations in Fig. [7j For the Type-I elongation, the 
increase in el energy and the decrease in hi energy im- 
plies a blue shift of E ff irrespective of the orientation of 
the elongation. This is same as earlier calculated for the 
flat QD. However, whereas for the flat QD, E g red shifts 
for both the Type-II and the Type-IIv elongations inde- 
pendent of their orientations, in this case of the tall QD, 
the changes in E s depend on the orientation of the elon- 
gation. For the [110] elongations, E s red shifts, whereas 
it blue shifts for the [110] oriented elongations. 
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FIG. 10: The polarization dependent optical transition modes TE[ 110 ], TE[ I10 ], and TMpoi] are drawn as a function of the (a) 
Type-I, (b) Type-IIv, and (c) Type-II QD elongations. The figures are plotted using the same scales to facilitate easy mutual 
comparison, (d) The plots of the in-plane polarization anisotropy (Pol||) as defined by Eq.[2]are shown for the three type of the 
elongations, exhibiting an inverse quadratic dependence on rj, consistent with the Sheng et al? 1 i: . Four cases are also marked by 
using green ovals indicating that the two different types of elongations, with roughly similar values of Tj, exhibit similar values 
of the Pol||. 



2. Polarization properties of the tall QD 

The polarization properties of the tall QD are signifi- 
cantly different from the flat QD of the previous subsec- 
tion because of the two major differences in the hole wave 
functions even for the circular-base case, as evident from 
the comparison of Figs. [3] and [HJ all of the hole wave 
functions for the tall QD are oriented along the [110] di- 
rection and are confined inside HH pockets at the QD 
interfaces in contrast to the flat QD where the hole wave 
functions are either nearly symmetric at the QD center 
or [110] oriented. Thus, for the circular-base tall QD, 
we calculate TE[ 110 ] < TE [TlQ j, and overall TE and TM 
modes have smaller magnitudes compared to the flat QD 
due to the relatively smaller spatial overlaps between the 
electron and the hole wave functions. 

We also want to specify that for the tall QD, the 
lower lying hole wave functions have larger contribution 
in the ground state optical transition intensity when com- 
pared to the flat QDs where mainly el-hl transition is 
dominant^. The dependence on the elongation factor 
is also quite different with all the polarization modes 



TE, 



no], 



TE [110]> 



and TMjq 01 ] increasing irrespective of 
the type and the orientation of the elongation. The large 
increase in the TEr-^i mode for the elongation along the 
[n] is same as also calculated for the flat QD and is 
explained in the section IV-A-2. However, the small in- 
crease in the other TE mode is attributed to an increased 
electron-hole spatial overlap because of elliptical shape of 
the QD. 

Table [TT| provides the ratio of the LH components in 
the top five valence band states for the elliptical QDs 
with respect to the circular-base QD. In each case, we 
also provide the corresponding ratio of the spatial overlap 
between the electron and hole wave functions (elhi) along 
the [001] -direction as defined earlier for the flat QD. The 



increase in the TMr 00 i] mode is directly related to an 
increase in the LH mixing of the valence band states. 
However for a giving LH mixing the magnitude of the 
TM[ 001 ] is proportional to the spatial overlap between the 
electron and hole wave functions along the [001] direction. 

As the shape of the QD is elongated, the electron wave 
function el also gets elongated along the major axis of 
the ellipse as evident from Fig. [8| This increased spread 
of el wave function results in a relative increase in its 
spatial overlap with the hole wave functions for the el- 
liptical QDs with respect to the circular-base QD as also 
noticeable from the values of elhi provided in Table |TT| 
For rj > 1.0 ([110] elongation), the LH mixing increases 
and along with an increase in elhi is responsible for the 
enhanced TMrpoi] mode. The [110] elongation (rj < 1.0) 
slightly reduces the LH character of the valence band 
states, however as the hole wave functions are oriented 
along the minor-axis, so a large increase in the electron- 
hole wave function spatial overlap overcomes a small (< 
15%) decrease in the LH component and therefore causes 
a net increase in the TM^oi] mode strength. It should 
also be noted that for the same amount of elongation, 
rj > 1.0 results in larger TM[ 001 ] mode when compared 
to rj < 1.0. This is because of the fact that for the r/ > 
1.0, both the LH mixing and the spatial overlap increase, 
whereas for the 77 < 1.0 the LH mixing decreases and 
only the spatial overlap increases. 

In-plane polarization anisotropy: As shown in 
Fig. 10 'd), the dependence of Pol|| on the elongation fac- 



tor (rj) for the tall QD is also quite similar to the case 
of the flat QD. The calculated values of Pol|| for all of 
the three types of the elongations again follow inverse 
quadratic dependence. Furthermore, for the fixed height 
of the QD, similar values of rj exhibit nearly same in-plane 
anisotropy. We highlight four such cases in Fig. [5jd) us- 
ing green dotted circles where nearly same values of the 
Pol|| are calculated for the following sets of the values of 
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TABLE II: The ratio of the LH components in the highest five valence band states (hi, h2, h3, h4, h5) of the elliptical-shaped 
tall QD with respect to the LH components of the corresponding valence band states of the circular-base (r/) QD for a few 
selected values of the elongation factor, rj. For each case, the corresponding ratio of the electron-hole wave function spatial 
overlap along the z = [001] direction given by elhi = |(i/; e i|z|-!/)hi}| is also provided. 
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FIG. 11: The plots of the DOP[no] an d DOP[y 10 j are shown as a function of the (a) Type-I, (b) Type-IIv, and (c) Type-II 
elongations. A large decrease in the value of the DOP[,i] along the minor-axis of the elliptical QDs is calculated. 



77: (1) 0.64 in Type-IIv and 0.67 in Type-II, (2) 0.82 in 
Type-II and 0.8 in Type-I, (3) 1.22 in Type-II and 1.25 
in Type-I, and (4) 1.565 in Type-IIv and 1.50 in Type-II. 
Therefore, we extend our conclusion of the previous sub- 
section that the in-plane polarization anisotropy (Pol||) 
only depends on the value of the 77, irrespective of the 
type and orientation of the elongation for both fiat and 
tall QDs, provided the height of the QDs is kept fixed. 

Tuning of DOP\^y. To conclude this discussion 
about the polarization response of the tall QD, we com- 
pare the values of the DOP[ 110 ] and DOPjy 10 ] for the 

[a), (b), and (c). 
for the circular-base (rj = 1.0), 



three types of elongations in Figs. 
Since TE [110] > TE [Tl0] 
so a much larger difference («13%) is present between the 
DOP[no] an d DOP[y 10 j for r\ — 1.0, as compared to only 
«0.2% difference for the flat QD. This difference is in 
agreement with the previous comparisorP^ between the 
similar flat and tall QDs and is attributed to the stronger 
orientation and confinements of the hole wave functions 
for the tall QDs. 

The elliptical shape of the tall QD once again reduces 
the value of the DOPp^j along its minor-axis. This is 



mainly because the TM[ oi]-mode increases for the elon- 
gated QDs whereas the corresponding TE-mode along 
the minor-axis does not increase much. In contrast to the 
case of the flat QD, the elliptical shape, however, does 
not reduce the value of the DOPp^j along the major-axis 
which in fact first slightly increases and then remains 
nearly constant. This is because for the flat QD, the TE- 
mode is much larger than the TM-mode and hence the 
values of the DOP^j are not very sensitive to changes 
in the TM-mode. However, for the tall QD, the TE- 
modes are comparatively smaller in magnitude due to 
the smaller electron-hole spatial overlaps, so the values 
of the DOPr-^i become more sensitive to the changes in 
the TE- and TM-modes. For the small values of r/, a 
larger increase in the TE-mode produces slight increase 
in the values of the DOP^j. For the large values of 77, 
the increase in the TM-mode also becomes important and 
hence the values of the DOP^j do not show any further 
increase. 

To summarize this subsection, we find that overall a 
much larger tuning of the polarization response is pos- 
sible by elongating the tall QDs. The largest reduction 



16 



(«51%) in the value of the DOPr-^i is calculated for the 
Type-II elongation at rj = 0.54. Since the red shift of 
the optical wavelength and the isotropic polarization re- 
sponse, both, are desired for the design of the optical de- 
vices operating at telecommunication wavelengths (1.3- 
1.5 /an), our model calculations suggest that the Type-II 
[110] elongations are more suitable as they fulfil both re- 
quirements (see Figs.[7^c), (f) and Fig. [TTJc)). It should 
also be noted that the confinements of the hole wave func- 
tions at the interfaces for the circular-base tall QDs re- 
duce the oscillator strengths by an order of magnitude 
when compared to the circular-base flat QDs. However, 
with the elliptical shapes, the oscillator strengths increase 
and even become comparable to the flat QDs, in partic- 
ular for the large Type-II elongations. 



C. Vertical Stack of Nine QDs (9-VQDS) 

Vertical stacks of QDs (VSQDs) have shown great po- 
tential for t uning of the polarization properties. Re cent 
experiment j 6 ! 9 ! 29 !^ and theoretical investigations^' have 
demonstrated that an isotropic polarization response can 
be realized by geometrical engineering of the VSQDs. In 
this subsection, we study a vertical stack of closely spaced 
nine QD layers (9- VSQDs) as shown by the schematic di- 
agram of Fig. [IJcVThis 9- VSQDs has been a topic of 
the recent studies 312121 due to its significant technological 
relevance to achieve isotropic polarization for the imple- 
mentation of the QD based SOA's. The optimized geo- 
metrical parameters of the 9- VSQDs are chosen directly 
from the experimentP, so that our results remain relevant 
to the experimental community. 

We have recently showrP by experimental PL mea- 
surements and theoretical calculations that the 9- VSQDs 
can exhibit TM[ 00 i] > TE[ 110 j leading to DOP[ 110 ] < 0. 
However, a significant anisotropy in the in-plane TE- 
mode was measured resulting in TEp- 10 ] > TMpoi] and 
DOPjj 10 ] > 0. Similar anisotropies in the DOPp 



were independently measured by alonso- Alvarez et o/J*«l 
and Humlicek et aiP^, which they were unable to ex- 
plain. Our multi-million atom simulations^ assuming 
circular-base for the QDs qualitatively explained that 
this anisotropy (DOPjuo] DOP[y 10 j ) is due to a strong 
confinement of the hole wave functions at the interfaces 
of QDs (similar to the case of the tall QD of the pre- 
vious subsection) which tend to align along the [110]- 
direction, and thus significantly reduce the TE[ 110 ]-mode. 
The TE[Y 10 ]-mode, on the other hand, does not observe 
any such decrease. The small increase in the TMfooi] 
mode due to the relaxation of the biaxial str ain, in par- 
ticular around the center of the 9-VSQDs^, is also not 
sufficient to overcome the TEry 10 j mode and thus the 
DOPjy 10 ] remains considerably larger than zero. 

The good agreement of our theoretical results with the 
experimental PL measurements even for an ideal circular- 
base 9-VSQDs shape leads to a fundamental question 
that how much is the contribution from the realistic 



shapes which are normally elongated. Here we system- 
atically elongate the QD layers inside the 9-VSQDs base 
and analyse its impact on the electronic and polarization 
properties. In contrast to the single QDs of the previous 
two subsections, where the QD base elongations result in 
a significant tuning of the DOPp^j, the magnitude of the 
DOPj^i for the 9-VSQDs is relatively insensitive to the 
value of rj. However, the sign of the DOPp^i is a strong 
function of the orientation of the elongation, and even a 
very small elongation (0.5-1.0 nm) is sufficient to control 
the sign of the DOPr-^i. Furthermore, we explore the 
possibility to achieve DOPr-^i < for both [n] = [110] 
and [TlO] by elongating the 9-VSQDs along the [110] di- 
rection. Our calculations predict that such a scenario is 
not possible due to a very high sensitivity of the DOPr^i 
with the value of r), and therefore the value of the DOPp^i 
may be reduced below zero for only one of the two spatial 
directions. 



1. Electronic properties of the 9-VSQDs 

Lowest three conduction band energies: 
Fig. [l2|a) plots the lowest three conduction band 
energies (el, e2, e3) as a function of the elongation 
factor (rj) for the Type-II elongation. The corresponding 
wave functions are also shown in Fig. |13| indicating 
that all of the lowest three conduction band states have 
s-type symmetry. Due to the strong coupling between 
the closely spaced QD layers, these electron states are 
hybridized over multiple QD layers^. Only a very small 
increase (less than 10 meV) in the energies of el, e2, and 
e3 is calculated as the 9-VSQDs is elongated. Due to 
the s-type symmetry, the effect of the Ad is negligible. 
We also find that in contrast to the single QDs, the 
energies of these molecular electron states are insensitive 
to the AV, as also confirmed by Fig. |l2"|c) where no 
shift in the el, e2, and e3 energies is calculated when 
the impact of strain is excluded. This is due to the fact 
that the electron energies are spread over multiple QD 
layers as well as occupies the GaAs spacers in between 
them, so a small decrease in the QD volume due to the 
Type-II elongation only negligibly impact their energies. 
Our calculations find that the hydrostatic component 
of the strain (e xx + e yy + e zz ) slightly increases for the 
elongated 9-VSQDs and results in a small upward shift 
in the conduction band energies. Also, the shifts in the 
energies are nearly equal for el, e2, and e3, and therefore 
the separation between them remains nearly unchanged. 

Highest five valence band energies: Fig. [l2^b) 
plots the energies for the highest five valence band states 
(hi, h2, h3, h4, and h5) as a function of the Type-II 
elongation. All of the hole energies increase for both, 
the [110] and the [TlO], elongations. 



Fig 



13 



shows the 



top view of the hole wave functions for 77=0. 57, 1.0, and 
1.75. For the circular-base 9-VSQDs, all the hole wave 
functions are aligned along the [110] direction due to the 
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FIG. 12: The plots of (a) the lowest three conduction band energies (el, e2, e3) and (b) the highest five valence band energies 
(hi, h2, h3, h4, h5) for the 9-VSQDs as a function of the Type-II elongation factor, rj. All of the energies increase irrespective 
of the orientation of the elongation, (c, d) The plots of the electron and hole energies as in (a) and (b), but without including 
the effect of strain. In the absence of strain, only Ad and AV contribute in the energy shifts, (e) The plots of the biaxial 
strain component (e xx + e yy — 2e zz ) along the [110] direction through the center of the 9-VSQDs with circular-base (77=1) and 
the [110] elongated base (r/=l. 75). (f) The plots of the biaxial strain component (e xx + e yy — 2e zz ) along the [110] direction 
through the center of the 9-VSQDs with circular-base (T]=l) and the [110] elongated base (?7=0.57). 
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FIG. 13: Top views of the wave function plots for the lowest three conduction band (el, e2, and e3) and the highest five valence 
band (hi, h2, h3, h4, and h5) states are shown for the circular-base 9-VSQDs and the two selected Type-II elongations of 
the 9-VSQD. The intensity of the colors in the plots represent the magnitude of the wave functions, with the dark color color 
indicating the smallest magnitude and the light green color indicating the largest magnitude. The boundaries of the QDs are 
also shown to guide the eye. 



presence of the HH pockets, simil ar to the case of the 
single QDs with large aspect ratioi^^l anc [ th e bilayer 
QDsP 

For the elliptical 9-VSQDs, the hole wave functions 
align along the major- axis. This is due to the larger bi- 



axial strain along these directions which pushes the HH 
pockets towards higher energies and hence the hole wave 
functions residing in these pockets also move towards 
higher energies. Fig. 12 'e) and (f) plots the biaxial strain 
components (e xx + e yy — 2e 22 ) through the center of the 
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FIG. 14: (a) The plots of the polarization dependent TE and TM modes as a function of the Type-II elongation factor 77 for 
the 9-VSQDs. (b) Plot of the P01|| as a function of the Type-II elongation factor for the 9-VSQDs. The plots indicate a high 
degree of in-plane anisotropy for the 9-VSQDs, which is very sensitive to the elongation factor, (c) Plots of the DOPjuoi and 
DOP[y 10 j as a function of the Type-II elongation factor for the 9-VSQDs. For the circular-base 9-VSQDs, 77=1.0, the DOP[no] 
is negative and the DOP[j 10 ] is positive. As we [110]-elongate the 9-VSQDs by 1 nm (marked by blue oval), the DOP[y 10 j 
becomes close to zero, however the DOP[no] drastically increases to +0.60, suggesting that it is not possible to simultaneously 
engineer both DOPs below or close to zero for this 9-VSQDs using [110] elongation engineering. 



9-VSQDs along the [110] and the [110] directions, respec- 
tively. The large increase in the biaxial strain is clearly 
evident for the elliptical 9-VSQDs when compared to the 
circular-base case. Since the hole energies are pushed up 
by an increase in the biaxial strain component, this shift 
dominates the small downward shift due to a small in- 



2. Polarization properties of the 9- VSQDs 



crease in the hydrostatic component (e a 



a result, all of the hole energies are shifted towards the 
higher values as evident from the Fig. [T2^b) . The contri- 
butions from Ad and AV are also very small, which is 
confirmed from Fig. 12 'd), where the strain effect is ex- 
cluded, and the hole energies move towards lower values 
due to the combined shift induced by Ad and AV. 



In Fig. |l2[b) , a relatively smaller increase in the hole 
energies for the [110] elongation as compared to the [110] 
elongations is due to the fact that all of the hole wave 
functions are initially oriented along the [110] direction 
for the circular-base (77 = 1.0), and they go through a 
90° rotation to align along the [110] direction for the 
[110] elongation. It should also be noted that the hole 
energy separations reduce as a function of the elongation 
factor (77) for the 9-VSQDs, suggesting even enhanced 
contributions from the lower lying valence band states in 
the ground state optical intensity measured at the room 
temperature for the elliptical 9-VSQDs. 



Optical gap energy, E g : From Figs. 12 a) and (b), 



both the electron and hole energies increase as a function 
of the elongation, however the hole energies have larger- 
slope and therefore the optical gap energy, E ff = ei - hi, 
decreases as a function of the elongation factor (77). We 
calculate a maximum red shift of «30 meV in E g for 77 
= 0.57. 



Fig. 



TE [110] 



14 



a) plots the polarization dependent TE[ 110 ], 
and TM[ooi]-modes as a function of the Type-II 
elongation. For the circular-base (77=1.0) case, TE[ 110 ] < 
TMjooi] and TE [I10] > TEr U0 ]. 

When the base of the 9-VSQDs is elongated, the in- 
crease in the biaxial strain component (see Figs. |12[ e) 
and (f)) increases the splitting between the LH and the 
HH bands, which will decrease the TM[ oi]-mode. How- 
ever, the small increase in the TM[ 001 ]-mode as plotted in 
the Fig. [l4]^ a) , is due to the hole wave functions confin- 
ing towards the middle of the 9-VSQDs for the elliptical 
shapes where the HH/LH intermixing is larger as com- 
pared to its edgeiPl For example, by comparing the 77= 1.0 
and 77=0.57 cases, we find following shifts in the spatial 
positions of the top-most five hole wave functions within 
the 9-VSQDs: hi moves from the QD layer 2 to the QD 
layer 3, h2 moves from the QD layer 3 to the QD layer 
5, h3 moves from the QD layer 2 to the QD layer 4, h4 
stays in the QD layer 8, and h5 moves from the QD layer 
3 to the QD layer 5; here the QD layers in the 9-VSQDs 
are numbered from 1 to 9 starting from the bottom to- 
wards the top, as mentioned in the schematic diagram 
of Fig. [ljc). Similar trends in the spatial confinements 
of the hole wave functions are observed for other values 
of 77 which are responsible for the small increase in the 
TM[ooi]-mode as a function of 77. 

The magnitude of the TE^i modes is calculated to be 
very sensitive to the elongation of the 9-VSQDs. Even for 
a very small [110] elongations (< 1 nm), the TErno] mode 
quickly increases above the TM[ 00 i]-mode. For example, 
for 77=1.095, the TEjY 10 j reduces to become close to the 
TMrooii * but the TE[ 110 ] has already increased by a factor 
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FIG. 15: Normalized Polar plots for the 9-VSQDs system with (a) circular-base, (b) Type-II [110]-elongation (77=1. 095), and 
(c) Type-I [110]-elongation (?7=1. 047). In each case, the polarization direction of the incident light (n) is kept along the [In- 
direction (red squares) for the TEj Il0 j-mode and along the [110]-direction (blue circles) for the TE[ 110 j-mode. The polar plots 
based on the cumulative sum of the optical transition strengths between the lowest conduction band state and the highest five 
valence band states are drawn with respect to the angle 9 between the [001] -direction and either the [110] or the [110] directions. 



of «13. 

It should also be noted that while the elliptical shape 
increases the TEn^-j mode along its major-axis, it has only 
very little impact on the TE-mode along its minor-axis. 
When the 9-VSQDs is elongated along the [110] direc- 
tion, the TE[Y 10 ] increases, but the TE[ 110 j mode remains 
nearly unchanged. Similarly for the [110] elongation, the 
TErY 10 j-mode quickly decreases and then remains nearly 
unchanged for 77 > 1.2. 

In-plane polarization anisotropy: Fig. [T4^b) plots 
the in-plane polarization (Pol| | ) defined by Eq. [2] as a 
function of the Type-II elongation factor (77). The large 
magnitudes of the Pol|| for both [110] and [110] elon- 
gations suggest a high degree of the in-plane polariza- 
tion anisotropy for the 9-VSQDs. Even for the per- 
fectly circular-base (77=1. 0), Pol|| is ~0.82, indicating 
that TEjY 10 ] » TE[ 110 j. Any elongation along the [110]- 
direction further increases this anisotropy. The [110]- 
elongation sharply increases TE[ 110 ]-mode and changes 
the sign of Pol||. This is because of the 90° rotation 
of the hole wave functions (see Fig. 13 1. Even an elon- 
gation as small as of 1 nm along the [110]-direction 
can change the value of the Pol|| from +0.82 to -0.51. 
Therefore, we conclude that the 9-VSQDs exhibits highly 
anisotropic in-plane polarizations. Similar in-plane po- 
larization anisotropics were measured by Alonso- Alvarez 
et a/P^ and Humlicek et alW^ for the vertical QD stacks. 

It should also be noted that whereas the in-plane po- 
larization (Pol||) for the single QDs, irrespective of their 
AR, exhibits an inverse quadratic relation with respect 
to 77 (see Figs.^d) and [ToJ^d)) , it demonstrates nearly a 
step function like dependence on 77 for the strongly cou- 
pled 9-VSQDs. 

Tuning of DOP^: Fig. W^c) plots the DOP [110] 

Dime 



For the ideal circular-base case (77 = 1.0), DOP[no] 
and DOP[j 10 j have values of -0.45 and +0.6. This 
is in qualitative agreement with the experimental PL 
measurements^ and leads to a question that how much 
impact would be from a realistic elliptical shape. Our 
calculations show that the orientation of the base elon- 
gation determines the sign of the DOP^j, whereas the 
magnitude of the elongation (value of 77) has a very little 
impact on the magnitude of the DOPp^j. This is clearly 
evident from Fig. 14 [c) for 77 < 1.0 and for 77 > 1.2, where 



and DOP[Y 10 j as a function o: 



Type-II elongation. 



a very small change in the magnitude of the DOPp^j is 
observed as the value of 77 is changed. 

The strong dependence of the sign of the DOPp^j on 
the orientation of the elongation is highlighted by using 
an oval in the Fig. 14 [c), where even for a 1 nm [110]- 
elongation, the DOP[ 110 ] drastically changes its sign from 
-0.45 to +0.6. This large change in the value of DOPp^j 
for 1.0 < 77 < 1.2 is remarkable as it indicated that only a 
very small shape asymmetry is capable of overcoming the 
effect of atomistic symmtery lowering effect. This also 
implies that the elliptical shape of the 9-VSQDs can not 
be exploited to simultaneously engineer both DOP[no] 
and DOPrY 10 i below zero. 

We want to highlight that although a tuning of the 
DOPp^j over a wide range of values is possible by the 
elongation of the single QDs, it remains relatively insensi- 
tive with respect to the magnitude of 77 for the 9-VSQDs. 
Therefore, we expect that the elongation of the 9-VSQDs 
would not offer much improvement in its polarization re- 
sponse. 

Polar plots: The strong impact of the [110]- 
elongations on the polarization properties of the 9- 
VSQDs is further confirmed in Fig.[l5]by comparing the 
normalized polar plots for (a) circular-base (77=1. 0), (b) 
Type-II 1 nm [110]-elongation (77= 1.095), and (c) Type-I 
1 nm [110]-elongation (77=1. 047). In each case, two polar 
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plots are drawn: (i) as a function of the angle 9 between 
the [001]-direction and the [110]-direction for the TEr 110 ] 
(blue circles); (ii) as a function of the angle 9 between 
the [001]-direction and the [110]-direction for the TErY 10 , 
(red squares). 

As the 9-VSQDs is elongated, a 90° rotation of the 
polar plots is calculated for the TEr 110 ]-niodes in both 
the Type-I and the Type-II cases. This clearly suggests 
that both [110] elongations will result in TEj 110 ]-mode > 
TMpoij-mode. For TE[ Tl0 j-mode, the polar plot rotates 
anti-clockwise by «30° and w45° for the Type-II and the 
Type-I elongations, respectively. This causes a reduction 
between the relative magnitudes of the TErj 10 i-mode and 
the TM[ooi]-modes, thus reducing DOPp; 10 ] from +0.6 to 
0.102 in (b) and to -0.07 in (c). 

Geometry of the experimentally grown 9- 
VSQDs: The above discussion about the strong de- 
pendence of the polarization properties on the [110]- 
elongations allows us to theoretically probe the geomet- 
rical shape of the 9-VSQDs as grown by Inoue et alW. 
It was reported that the 9-VSQDs are not isotropic and 
the TEM images suggested very little anisotropy in the 
lateral extentp2, possibly a [110]-elongatiorP2l. Our multi- 
million-atom calculations show that the polarization re- 
sponse is very sensitive to the elongation factor (rf) and 
even a 0.5-1.0 nm [110]-elongation increases DOP[ 110 ] 
above zero. Therefore according to our model results, the 
experimentally measured DOP[ 110 ] = -0.6 implies that 
the shape of the 9-VSQDs can only have [110] elongated 
base which confirms the findings from the TEM images^!. 
It should also be noted that as the 9-VSQDs studied 
here has pure InAs QD layers, so our finding does not 
contradict with the conclusions of Mlinar et a/P3 where 
they report that for the alloyed InGaAs QDs, the alloy 
random configurations may significantly impact the po- 
larization properties and make the correlation between 
the measured polarization response and the QD geome- 
try unreliable. 

V. SUMMARY AND CONCLUSIONS 

We have performed multi-million-atom simulations to 
understand the impact of the elliptical shapes on the elec- 
tronic and polarization properties of the single and the 
multi-layer vertical stacks of InAs QDs. The comparison 
between a flat QD and a tall QD, having aspect ratios 
of 0.225 and 0.40 respectively, reveals drastically differ- 
ent electronic and polarization properties as a function 
of their base elongation. The key outcomes of the com- 
parison are: 

(i) The quadratic component of the piezoelectric poten- 

tial completely cancel the linear component inside 
the tall QD region, whereas only partial cancella- 
tion occurs for the flat QD. 

(ii) Although the stain and the piezoelectric potentials 

are drastically different for the fiat and tall QDs, 



the lower electron p-state (e2) is oriented along the 
[110] direction for both systems. 

(iii) The hole wave functions are confined inside the flat 
QD mainly at its center, whereas they are con- 
fined at the QD interfaces inside the HH pockets 
for the tall QD. This leads to a reduction of the 
oscillator strengths for the circular-base tall QD 
by approximately an order of the magnitude due 
to smaller electron-hole wave function spatial over- 
laps. For the elliptical-shaped tall QDs, the oscilla- 
tor strengths increase and even become comparable 
to the flat QD for some values of the elongation fac- 
tor (ij). 

(iv) The Type-I elongation, irrespective of its orienta- 
tion, blue shifts the optical gap energy (E 9 ) for 
both, the flat and the tall QDs. The Type-II and 
Type-IIv elongations red shift E g for the flat QD 
irrespective of the elongation direction, whereas for 
the tall QD the shift in E fl strongly depends on the 
orientation of the elongation: red shift of E s for the 
[110] elongation and blue shift for the [110] elonga- 
tion. 

(v) The elliptical shape of the flat QD always improves 

its polarization properties by reducing the value of 
the DOPpf], whereas it only reduces the DOPp^j 
along the minor-axis for the tall QD. 

(vi) The elliptical shape of the tall QD allows a tun- 
ing of the DOPp^-j over a much wider range, when 
compared to the flat QD. This property can be fur- 
ther exploited in large stacks of strongly coupled 
QD layers, where essentially very high values of the 
ARs can be achieved. 

Although the understanding of the single layers of the 
QDs provides significant physical insight of the impact of 
the shape asymmetry, they do not lead to isotropic polar- 
ization response (DOPp^j ~ 0) for the QD base elonga- 
tions of up to 6 nm studied in this paper. Therefore, we 
prbextend our study of the elliptical shapes to the exper- 
imentally reported vertical stack of nine QDs (9-VSQDs) 
which has demonstrated DOP[ 110 ] < 0. The key features 
of our analysis about the elliptical 9-VSQDs are: 

(i) In contrast to the single QD layers where the ellipti- 

cal shape only very slightly reduces the magnitude 
of the biaxial stain, a significant increase in the 
magnitude of the biaxial strain is calculated for the 
9-VSQDs. Therefore, the shifts in the hole energies 
are dominated by the changes in the biaxial strain, 
rather than AV and Ad. 

(ii) The hole wave functions are confined inside the HH 

pockets at the interfaces of the 9-VSQDs, similar 
to the case of the tall QD. This introduces a large 
in-plane polarization anisotropy. 
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(iii) While the elongations of the single QDs result in the 
tuning of DOPp^i over a wide range, the magnitude 
of DOPp^j is largely insensitive to the magnitude 
of ij for the 9-VSQDs. Therefore, we conclude that 
the elliptical shapes of the 9-VSQDs do not provide 
any noticeable improvement in the polarization re- 
sponse. This is clearly evident from Fig. 14 'c) for 
T] < 1.0 and for r) > 1.2. 

(iv) Our calculations show that the sign of DOPp^-j is 
very sensitive to the orientation of the elongation. 
Even a very small variation of r\ from 1.0 is capable 
of controlling the sign of the DOPp^j. We find that 
the in-plane polarization (Pol||) roughly follows a 
step function like abrupt dependence on 77, as com- 
pared to an inverse quadratic dependence for the 
single QDs. Such a large in-plane anisotropy of the 
polarization allows to accurately predict the shape 
elongation of the 9-VSQDs studied in this paper, 
in agreement with the TEM findings. 



In summary, we have presented a detailed analysis of 
the dependence of the polarization properties as a func- 
tion of the elongation factor 77 that would serve as a guid- 
ance to engineer the geometry parameters for the tuning 
of DOPpyf] from semiconductor QDs, which is a critical 
design parameter for several challenging applications. 
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